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The pseudogap phase of the cuprate superconductors is argued to be character- 
ized by a hidden broken symmetry of d-wave character in the particle-hole channel 
that leads to staggered orbital magnetism. This proposal has many striking phc- 
nomenological consequences, but the most direct signature of this order should 
be visible in the neutron scattering experiments. The theoretical underpinning of 
these experiments is discussed. 



1 Introduction 

It has been proposed that the pseudogap state is a phase with a broken sym= 
metry and an order parameter ,Em not a state with a merely fluctuating order u 
The distinction between the two is enormous. The order parameter, which is 
due to a particle- hole condensate with internal "angular momentum" 2, will 
be termed the d-density wave (DDW). We have predicted that the present 
day neutron scattering experiments have enough resolution to provide a di- 
rect evidence of this order. And, indeed, there are some preliminary neutron 
scattering measurements in the bilayer YBa2Cu306.|Sr-t na t appear to have 
observed the theoretically predicted order parameter .HeI Thepe are other ex- 
periments on YBa2Cu30e.5, which observe a magnetic signaljs but appear to 
be from remnant spins left over from the undoped antiferromagnet. While the 
experimental situation must be resolved, we shall focus merely on the theoret- 
ical predictions because a precise analysis is necessary to make sure that what 
is being observed is not an artefact. The fundamental signature of this order 
parameter is an elastic Bragg peak centered at two-dimensional wave vector 
Q = (7r/a, 7r/a), where a is the lattice spacing of the CuO-lattice. Since the 
magnitude of the staggred orbital fields is of order 10 G, the experiments are 
difficult. _ 

The neutron scattering from DDW was addressed previouslyQ but the 
polarized neutron scattering and the important question about the current 
form factors were not. The current form factors play an important role in 
polarized neutron measurements, without them one may be led to physically 
incorrect results. Moreover, the crucial connection to the pseudogap phase was 
not made, nor was it made clear that such a phase truly survived beyond the 
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large- N mean field theory due to presumed gauge fluctuations. There is also 
a proposaB of a circulating current phase, which does not break translational 
symmetry. The similarity with the DDW state is merely superficial. At the 
very least, the signature, on robust symmetry grounds, is a neutron signal at 
Q = which is fundamentally different. We will have nothing to say about 
this proposal, as it is not germane to the present set of experiments. 

2 Order parameter 

The singlet DDW is defined by the order parameter $q in terms of a particle- 
hole condensateQ 

{c af (k + Q, i)c p (k, i)> = (cos k x a - cos k y a) 8" p , (1) 

where c, p are spin indices. The order parameter breaks time reversal, parity, 
translation by a lattice spacing, and the rotation by tt/2, although the product 
of any two of these is preserved. Interestingly, DDW does not modulate charge 
at all, but it modulates currents. The reason it is still called a density wave 
is because it is a particle-hole condensate, and d-wave refers to the internal 
form factor of the particle and the hole, which is (cos k x a — cos k y a). We 
expect that the magnetic field generated by the circulating currents should 
be proportional to the pseudogap.0 A simple estimate results in a tiny field of 
order 10 G due to currents circulating in a CuO plaquette as shown below: 



Figure 1. Circulating current pattern in the model where the current is carried by 5-function 
wires. 



Although this simple model captures the correct symmetry, it is cruicially 
deficient in predicting the correct neutron scattering intensities, as we shall 
see. It is a priori clear that the flow of currents cannot be along infinitesimally 
thin wires, but it must spread out, both because of finite sized atomic orbitals 
and because of electron-electron interactions. One can model the more real- 
istic current distributions by the form factors a(q) and (3(q) in the Fourier 
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transform of the current distribution (j(q)) without destroying the symmetry 
of the order parameter (q • (j(q)) = 0), where 
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Here, qy = (q x ,q v ), /(q) = sin(2|-), d is the separation between the two 
CuO-planes within a bilayer complex, and G is a reciprocal lattice vector. 
The simplest choice is to take a(q) and /?(q) to be dependent on q z only, such 
that a(q z ) — @{q z ) ~ , as q z — > 0. The choice of this current density was 
discussed in Ref. 5. 



3 Experimental detection of DDW 

The DDW is usually hard to detect because 

^(c cr ^(k + Q)cp(k)) cx ^(cosfc x a — cosk y a) = (3) 

k k 

There is no net modulation of charge and spin that could be measured. Con- 
trast this with a d-wave superconductors (DSC) for which a similar momentum 
sum of the particle-particle condensate also vanishes. Nonetheless, because 
of broken gauge symmetry, Meissner effect follows, irrespective of the pairing 
channel. This is not possible for a DDW because the broken symmetry is 
in the Ising universality class, and experiments seeking to uncover such order 
must (a) be sensitive to spatial variation of kinetic energy or currents, (b) mea- 
sure higher order correlations of the charge or spin density as in 2-magnon 
Raman scalletring, nuclear quadrupole resonance, etcQ One nice feature of 
DDW is that a thermodynamic transition is possible in two dimensions due 
to the Ising nature of the order parameter. 

Impurities will mix d-wave order with the s-wave order proportional to 
the concentration of impuritiues. This may be an indirect way to reveal the 
hidden broken symmetry, because the corresponding s-wave order is none 
other than the conventional charge density wave (CDW). Another possibility 
is spin-orbit coupling, which will mix spin density wave (SDW) with DDW. In 
mean field theoryjlj there is some evidence that DDW may incommensurate 
for higher doping. If this is the case, the charge order must inevitably be 
mixed in, an idea that can be experimentally probed by scanning tunneling 
measurements inside the vortex core of the mixed DDW and DSC phase. The 
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detailed nature of the charge order appear to be highly nonuniversal, and 
certainly cannot be used as a diagnostic tool for theories of pseudogap. 



4 Neutron scattering 

4-1 Unpolarized neutron scattering 

For unpolarized neutron scattering, the differential scattering cross sections 
for Bragg scattering from an ordered array of orbital currents is given byrl 

where q now is the momentum transfer. We certainly do not know the form 
factors a(q) and /3(q) beyond their limiting forms, but a reasonable guess can 
be made to see the robust features of the intensity modulation as a function of 
the momentum transfer. For this purpose, we have chosen a(q) = f(q)e ™ , 
and /3(q) = f(q)e 51 ■ So, the intensity is parametrized by a orthorhom- 

_ r/ Qz_ \2 / qz \2i 

bicity parameter A(q) = e 11 w ; - 1 1 and /3(q). Note that all other factors, 
such as the magnitude of the order parameter, drop out if we normalize with 
respect to a reference Bragg-ceflection intensity. We then replace /3(q) by the 
well known Cu form factorJIJ which is likely to be an upper bound because 
the orbital currents are more spread out in real space than atomic orbitals. 
We choose qo — 2ir/d and qi = 0.275<7o- We shall see that the unpolarized 
intensity is not sensitive to the choice of the A-parameter. 

There is no question that the unpolarized neutron scattering intensities 
are dramatically diffrent for orbital currents in comparison to spins lying in 
the a-b plane as shown in Fig. |^. The scattering from spins pointing along 
the c-direction can be easily calculated and is similar to the scattering from 
orbital currents. These two cases can be distinguihed by going to higher order 
Bragg reflections such (H, H, L), for H > 1/2, for example (3/2, 3/2, L). The 
results for orbital currents fall off very rapidly in contrast to spins pointing in 
the c-direction. 

4-. 2 Polarized neutron scattering 

The Bragg scattering intensity for polarized neutrons from orbital currents is 
given by 

^^lai/^-qxCKq))! 2 , (5) 
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Figure 2. Normalized intensity, /, of Bragg reflections as a function of the perpendicular 
momentum transfer, q z . The shorthand notation is I(H,K,L) = I(2nH/a, 2irK/a,2iTL/c). 
The intensity from spins lying in the a-b plane is obtained from the form factor given in 
Ref. 11. 



where jl is the neutron spin, and \i) and |/) are the initial and the final states 
of the neutron. It is easy to see that if the scattering vector is parallel to the 
direction of polarization, the entire scattering is spin flip. This is useful to 
identify magnetic scattering from non-magnetic scattering and the intensity 
ratios should follow the same pattern as that of the unpolarized case. When 
the scattering vector q is perpendicular to the polarization, there can be 
both spin-flip and non-spin-flip scattering. This geometry contains additional 
information, which is the ratio of the spin-flip to non-spin-flip scattering. A 
common experimental set up is to polarize neutron along along the [1,1,0] 
direction as shown in Fig. [| so that the new x-axis is along [1, 1,0], and the 
new y-axis is along [0, 0, 1]. 

Then a simple calculation shows that the ratio of the non-spin-flip to 
spin-flip scattering intensity, NSF/SF, for Bragg reflection (H, H, L) is 



which vanishes identically if there is no orthorhombicity, that is, A(L) = 1. 
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Figure 3. The vertical field geometry. 



This is a nontrivial result for orbital current theory. 

For the choice of the parameters described earlier, corresponding to X(L — 
1) w 0.4, this ratio for the reflection (1/2, 1/2, 1) turns out to be w 1. Thus, 
while modeling of the orthorhombicity of the current distribution had very 
little effect on the unpolarizcd intensity, its effect on the polarized scattering 
is striking, turning the ratio (NSF/SF) to order 1 as compared to 0. 

5 Concluding remarks 

Other possible evidence against spin order is the existence of a spin gap due to 
the Ising universality class of the DDW order, and the small moment of order 
2 x 1Q~ 2 /ib corresponding to the field due to circulating currents. We made an 
educated guess for the microscopic nature of the current distribution. Clearly, 
we need a better understanding of it, although the ratio of the intensities 
as a function of the momentum transfer for unpolarized scattering is very 
insensitive to the modeling of the current distribution. Thus, we believe that 
our results are quite robust. 
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